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In the description of the interaction between electrons beyond the classical Hartree picture, bare 



^ ' exchange often yields a leading contribution. Here we discuss its effect on optical spectra of solids, 

♦~H^ comparing three different frameworks: time-dependent Hartree-Fock, a recently introduced com- 

■^^ . bined density-functional and Green's functions approach applied to the bare exchange self-energy, 

and time-dependent exact-exchange within time-dependent density-functional theory (TD-EXX). 

We show that these three approximations give rise to identical excitonic effects in solids; these 

effects are drastically overestimated for semiconductors. They are partially compensated by the 

usual overestimation of the quasiparticle band gap within Hartree-Fock. The physics that lacks in 

Jh I these approaches can be formulated as screening. We show that the introduction of screening in 

TD-EXX indeed leads to a formulation that is equivalent to previously proposed functionals derived 

from Many-Body Perturbation Theory. It can be simulated by reducing the long-range part of the 

Coulomb interaction: this produces absorption spectra of semiconductors in good agreement with 

experiment. 



INTRODUCTION 



Q ■ Density Functional Theory (DFT) [l| is one of the most widely used approaches for the calculation of material 

O \ properties that are determined by the electronic ground state. Since existing approximations for the in principle 

exact but in practice unknown exchange-correlation contribution are sometimes not sufficient to obtain the desired 

rvQ \ accuracy or even qualitative behavior, the search for better functionals is a continuous effort. The same is true 

I> concerning neutral electronic excitations that are accessible observables of time-dependent DFT (TD-DFT) [3|. The 

OO _ failure of most existing approximations for extended systems is in this case even more significant than in the case of 

^^ ■ the ground state Q- In both cases, it has been recognized that one may have to accept an additional complication of 

I ' the functionals in order to get reliable results. 
f~^ - One way to go are orbital-dependent functionals. A local Kohn-Sham (KS) [j] exchange-correlation (xc) potential 

VO i 'Vxc can be obtained from agiven non-local self-energy via the (linearized) Sham-Schliiter equation [5J . This optimized 
^D ' effective potential (OFF) 6] contains much of the important physics of the underlying self-energy. The most prominent 
approximation along this path is the so-called "exact exchange" (EXX) [TJ. The EXX potential wexx is obtained 
as the local counterpart for the non-local Fock operator, which is usually referred to as the "exchange" term. This 
potential is of course completely different from the original Fock operator: the OFF is local and constructed with KS 
i-rt ■ wavefunctions, whereas the exchange operator is non-local and constructed with Hartree-Fock (HF) wavefunctions. 
P5 It has been found empirically that EXX eigenvalue band gaps are closer to experimental quasiparticle band gaps 

O • than are local-density approximation (LDA) or HF ones. The agreement is very good for simple semiconductors 
[8|. However, an increasing underestimation for materials with wider band gap has been noticed '9l|. Although FXX 
contains important effects (it is devoid of self-interaction |3, as HF is), it is not meant to simulate the effects of 
Hartree-Fock. In particular, there is no simple interpretation for the KS eigenvalue band gap, no equivalent to the 
Koopmans theorem. Since KS band gaps are a priori not meant to reproduce experiment, the band gap agreement 
C^ ' may arise to a large extent from error cancellations (though still be useful). 

Instead, in the case of electronic excitations, the situation is different. Time-Dependent Hartree-Fock (TD-HF) 
can be understood as an approximation to Many-Body Perturbation Theory (MBPT), obtained from the latter by 
neglecting correlation. TD-FXX, on the other side, is a straightforward approximation of TD-DFT. Both MBPT and 
TD-DFT should in principle yield the same (correct) dynamical polarizability, that can be measured e.g. by optical 
absorption. 

HF and, more recently, EXX have been used quite frequently in calculations of real solids |l(i |ll|, sometimes 
augmented with approximate correlation functionals, like LDA [l3|. Much fewer exarnples exist instead for their 
time-dependent counterpart; TD-HF has been carried out for large band gap rnaterials Il3| . and, to our knowledge, 
only the absorption spectrum of silicon has been calculated within TD-EXX |lj| . Intuitively, one would of course not 



choose TD-HF to calculate, say, the absorption of silicon, since the strong screening in this material can be expected 
to drastically influence the electron-hole interaction. On the other hand, with EXX KS band gaps comparing so 
favorably to experiment with respect to their HF counterparts, one may hope for a similar improvement when going 
from TD-HF to TD-EXX. In fact, it was suggested in Ref. [ijthat the TD-EXX absorption spectrum of bulk silicon 
favorably compares to experiment. It is therefore worthwhile to elucidate the links between the two approaches and 
discuss differences in the various ingredients as well as in the expected and calculated results. This is the main aim 
of the present work. 

For the sake of clarity and completeness, we add to this comparison a third method, that is to some extent 
intermediate. This recently introduced approach, called p/G in the following, is in fact a combination of MBPT and 
the density- functional concept |1!tI |. As will be discussed below in the exchange-only approximation, it is situated 
on the same level as HF concerning the quasiparticle band gap, but close to TD-EXX concerning the electron-hole 
attraction. 

The three methods used in the present work are different for some aspects, but tightly linked for others. For a 
better understanding, we summarize in Tabled the differences and similarities of the three frameworks. The meaning 
of the notations and the acronyms will be made clear along this paper. 

Our mathematical comparison and numerical results (obtained for the example of bulk silicon) clearly show that 
in the case of solids all three methods are equivalent. This implies that TD-EXX in its present form is not suitable 
for the description of absorption spectra of semiconductors. We discuss hence the need for screening of the exchange 
interaction, and, more precisely, for the screening of the long-range components of it. By introducing the missing terms, 
we make the link with recently introduced and successfully used functionals derived from Many-Body Perturbation 
Theory [la, lla llll llM llS| j and we discuss which are the most crucially needed corrections. 

TABLE I: Schematic overview of the different approaches and corresponding observables, approximations and acronyms. 
Quantities are specified for the exchange-only case. 
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ELECTRON-HOLE INTERACTION IN THE QUASIPARTICLE AND IN THE DENSITY-FUNCTIONAL 

FRAMEWORK 



It is useful to first recapitulate the main features of the various approaches that will be compared. For the sake 
of clarity, we simplify the problem and do not distinguish between KS and HF single-particle wavefunctions in the 
following. In the numerical examples, we use LDA wavefunctions throughout, and furthermore, all the band structures 
(LDA, HF, EXX) differ solely by a rigid shift of the band gap. These assumptions are very reasonable for bulk silicon 
12(1. I2 ij (of course, the situation would be rather different in finite systems |2^). Atomic units are used throughout 
the present work. We often employ the widely spread many-body short-hand notation 1 = (ri, ii) and omit to specify 
spin explicitly. 

The neutral excitations of materials are described by the polarizability or density-density linear response function 
x(r,r',t — i'). This quantity expresses the linear response of the electronic density to variations of an external potential 
dU,^tir',t'y. 



X(r,r',i-i') = 



Sp{r,t) 



'5C/e.t(r',i') 



(1) 



Beside the response of independent particles Xo(r, r', t — t'), the full density-density response function further contains 
contributions stemming from the self-consistently induced potentials. Independent particle and interacting particle 
response functions are linked via a Dyson-like polarizability equation, symbolically: 



first of all is due to a self-consistently induced Hartree potential v'^^'^{r,t) 



(2) 



where the kernel k 

J dr'dt' [6vH{r, t)/6p{r' , t')] Sp{r' , t'). The induced Hartree potential contributes hence to the kernel k with the simple 

term [5vh{v, t)/Sp(r', t')] = v{r - r')6(t - t'). 

A similar variation has to be added for the exchange-correlation potential. In the TD-DFT framework, this leads 
to the exchange-correlation kernel fxci^.r' ,t — t') = [Svxc{Tc,t)/5p{r' ,t')], and Eq. (O is an integral equation for 
X{r,r',t-t'y. 



X(l,2) 



xr(l,2) 



Yd3d4xr(l,3)X 



TD-DFT 



(3,4)x(4,2). 



(3) 



Here, X (3,4) — w(3,4) + fxc{5,4). w(l,2) = t;(ri — r2)S{ti — ^2) stands for the instantaneous bare Coulomb 

interaction and Xq^ is the Kohn-Sham independent-particle response function. 

In the MBPT framework, the equivalent variation has to be calculated for the non-local exchange-correlation self- 
energy S. As a consequence, it is not possible to obtain straightforwardly a closed Dyson-like equation for the 
two-point density response function. This remains also true in the HF approximation, where S reduces to the - still 
non-local - Fock operator Y.^- Moreover, the exchange term has a simple dependence on the one-particle Green's 
function G'(rii,r2i^) [not so simple on the density p{r)]. Therefore, as in the general case, the derivation of a closed 
equation for the response function of the Hartree-Fock system leads to the introduction of four-point polarizabilities 
'^Xo and ^x- 



Q 



*X(1,2;1',2') 



^Xo(l, 2; 1', 2') + J dM4d5d6 4xo(l, 3; 1', 4)^/^(3, 6; 4, 5)^(5, 2; 6, 2') 



(4) 



[Contracting indices to '*x(l: 2; 1, 2) gives back the usual x(l, 2).] This equation is the Bethe-Salpeter equation J3L l23l| . 
For the Hartree-Fock case, the kernel '^K is [2S| 



'^K{3, 6; 4, 5) = 25(3, 4)5(5, 6)w(3, 5) - 5(3, 6)5(4, 5)w(3, 4). 



(5) 



The first v of the kernel stems from the variation of the Hartree potential (with a factor of 2 for singlet excitons), 
and the second v, from the variation [SYix/6G] of the Fock exchange operator S^;. ^XO: the 4-point independent- 
quasiparticle polarizability, is constructed here using HF eigenvalues, ^Xc?^ (see last column of Tabled. Note that 
beyond bare exchange, in most applications of the Bethe-Salpeter equation a statically screened Coulomb interaction 
W replaces v in the second term p| . This term gives rise to a direct screened electron- hole attraction instead of the 
unscreened one in the case of HF. 

Recently, Bruneval et al. [l5j introduced a general formulation that avoids the solution of the four-point Bethe- 
Salpeter equation. This approach (named p/G here) is based on the use of the density-functional concept within 



Hedin's equation of MBPT J24|. Of interest here, it yields an explicit formula for the two-point polarizability x that 
contains all exchange-correlation effects via the density-variation of the self-energy, STi/Sp, and via the explicit use of 
an independent-quasiparticle (QP) polarizability. The approach is between TD-DFT [because it leads to a two-point 
polarizability equation like Eq. Q] and MBPT (because the QP, instead of the KS, independent-particle response 
appears). Using the GW approximation for E J2J| and some further straightforward approximations, this approach 
leads to previously proposed equation for x that yields spectra in excellen t ag reement with the spectra calculated 
from Bethe-Salpeter equation, and consequently, with experimental spectra |l6l.ll7l ll8lll9| |. This equation is similar 
to Eq. ^, but recast into a two-point form: 



X(l, 2) = xo(l, 2) + y d3d4 xo(l, 3)J^''/«(3, 4)x(4, 2), (6) 

where the two-point kernel K^''^ reads 

KP"'{1,2) := «(1,2) + /(2)(1,2) (7) 

= «(l,2)-*yd3d4d5xo"'(13)G(34)G(53)^^. (8) 

For the present work, it is interesting to choose the Fock operator as an approximation for the self-energy in ST,/Sp, 
and consistently, to use the HF independent-particle response function Xo^ for xq in the polarizability equation 
(second column of Table P). The two-point kernel K''/'^ of the Hartree-Fock problem becomes, following Ref. llSl 

K''/«(l,2)^t;(l,2)+ /"d3(i4d5d6x"^"'(13)GHF(34)GHF(53)z;(45)GHF(46)GHF(65)xr"'(62), (9) 

where Ghf(1,2) is the HF Green's function and hence Y^x = iGnpv the Fock operator. Here, as in Ref. Ilfii we have 
used the approximation SG/Sp — —G{5G^^/5p)G c^ Gxq^G. 

One may wonder whether Eq. (jHJ using this approximate p/G kernel is indeed able to reproduce TD-HF. Fig. Q] 
compares the TD-HF (continuous curve) and p/G (dashed curve) calculated 12^ absorption spectra for bulk silicon. 
Of course, the results are far from any experiment: the HF direct band gap of silicon is 8.92 eV [see the independcnt- 
HF-QP result, dot-dashed curve obtained from Im(xo)], more than twice the experimental QP direct band gap of 
3.40 eV i2g;j. Also the electron-hole attraction is drastically overestimated due to the absence of screening, and a 
strongly bound exciton is formed inside the HF-QP band gap. Finally, QP and excitonic errors cancel to a large 
extent; the absorption spectrum falls in an energy region that is closer to the experimental one (circles) [25|, but the 
lineshape is of course completely wrong. This is to be expected and is not the point here. Instead, it is noteworthy 
to point out that the TD-HF [Eqs. Q and (0)] and p/G [Eqs. © and JHl] results are almost indistinguishable, as it 
was the case in previous findings when correlation beyond Hartree-Fock was taken into account [l3| and the effective 
interaction was therefore much weaker. This means that the kernel f^"^' in Eq. (0) simulates well the TD-HF hare 
electron-hole attraction. 

Let us now come to a fully TD-DFT formulation of the problem (see the first column of Table P) . This can also 
easily be written starting from the equations of Ref. [iSf indeed, that work showed that a differentiation of the time- 
dependent Sham-Schliiter condition 27] (that the TD-DFT and the MBPT time-dependent densities correspond) 
with respect to the density yields the TD-DFT polarizability equation, 

X(l, 2) = xr(l, 2) + I d2d3 xr(l, 3)i^™-DFT(3^ 4)^(4^ 2), (10) 

in which now the KS independent-particle response Xo"^ appears, and the corresponding kernel reads t^td-dft _ 
KP/'^ + {Xq Xo ^)- I'^ other words, the TD-DFT kernel has, with respect to the p/G one, an additional contribution 

/(^) := Xo "^ ~ Xq^ ■ If inserted in a Dyson-like polarizability equation, /'^^ transforms the KS independent-particle 
response Xo^ into the corresponding QP independent-particle response xo- It essentially opens the band gap from 
the KS to the QP one [Hil. 

When applied to the exchange-only case, one obtains hence k"^'^'^^"^ = w + Z^^-* + /'■^-' where /^^^ has the role to 
open the band gap to the HF one, and the bare electron-hole attraction /*^^) is the kernel of Eq. ^ [approximated 
e.g. by Eq. iPJ]. We call this approach the "non-linearized TD-EXX approach", as opposed to the standard TD-EXX 
approach that will be discussed in the next section. 

It is quite obvious to see that Eq. (|10|l and Eq. ® yield identical results, as we have also confirmed numerically 
(without displaying the result here). This leads to an important conclusion of this section, namely: non-linearized 



Figure 1, Bruneval et al., J. Chem. Phys. 
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FIG. 1: Imaginary part of the macroscopic dielectric function of bulk silicon. Continuous curve: TD-HF result; dashed line: 
pjG result (see text) from Ref. ilfit dot-dashed curve: independent-HF-QP result (see text); circles: experiment [2RJ|. 

TD-EXX reproduces TD-HF. The band gap difference between EXX and HF is cancelled in the optical spectrum by 
the contribution /^^^ to the kernel, whereas the effect of the electron-hole attractions STf^/SG and f'--^^ are extremely 
close. 



TIME-DEPENDENT HARTREE-FOCK AND TIME-DEPENDENT EXACT EXCHANGE 



Let us now make the link to what is usually called "TD-EXX" . The starting point is the linearized Sham-Schliitcr 
equation, ||27l 



i / dl WEXx(l)Xo (1,2) = / dldS / d4 Gexx(1, 3)I],(3, 4)Gexx(4, 2) 



(11) 



where only EXX KS quantities are used to build response functions and the Fock operator E^; — iGexxv (the solution 
of the static version of this equation would yield the static OEP potential ■yEXx)- The functional derivative with 
respect to the density of f exx has a contribution that stems from the derivative of J^x (this is explicitly shown in 
appendix J , 

/(2)'""(1, 2) - / d3dM5d6 Xo''''"'(l, 3)Gexx(3, 4)Gexx(5, 3)^(4, 5)Gexx(4, 6)Gexx(6, 5)xo''''"'(6, 2). (12) 



This term is very similar to f^"^' in Eq. @. Since we do not distinguish single-particle wavefunctions, the only 
difference lies in the eigenvalues used to build all the Green's functions, and xo- 



Figure 2, Bruneval et al., J. Chem. Phys. 
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FIG. 2: Imaginary part of the macroscopic dielectric function of bulk silicon. Continuous curve: pjG result (see text) using 
Z'^'; open circles: pjG result using /(^j.iin j^j^jg corresponds to TD-EXX (see text)]; dot-dashed curve: using jf^)'''" and the 
modified Coulomb interaction of Ref. VVM. 



Actually, it turns out that this difference is not significant. Fig.|21shows a comparison of the pjG spectrum calculated 
using /'^' and /'^''''", respectively (continuous line and open circles). The reason for this perfect agreement is a strong 
cancellation between energy denominators in the inverse response functions and the two Green's functions in terms 
of the type GGx^^ ■ 

It is now important to notice that /(2)>iin jg nothing else but the electron-hole attraction term of TD-EXX. This 
term corresponds precisely to the terms Il\ and ll\ in Ref. ^ (see appendix „ for a detailed derivation). Hence, wt 
find a strongly overbound exciton from the TD-EXX electron-hole attraction. 

The rest of the terms that stem from the derivative of Eq. Hll|l is the linearized version of f^^\ It reads (see 
appendix ^ 

/(^)'""(1, 2) = / d3d4d5d6 Xo''''"'(l, 3)Gexx(3, 6)Gexx(6, 4) [S,(4, 5) - 5{4, 5)wexx(4)] Gexx(5, 3)xo''''"'(6, 2) 
-f / dMMbdG Xo'^'^'^h 3)Gexx(6, 3)Gexx(3, 4) p,(4, 5) - (5(4, 5)t;Exx(4)] Gexx(5, 6)xo''''"'(6, 2). (13) 



As shown in appendix^ this expression (again assuming that KS and HE wavefunctions arc equal) corresponds to the 



terms iJ|- and ij| of Ref. ]2£ 



Z*^^' has the difficult task to shift the whole independent-particle spectrum above the HE band gap, and it turns 
out that it is more delicate to linearize this contribution than /'■^K It is clear that the linearization of /^^^ cannot, 
by miracle, cancel the overestimate of the exciton binding in /(^^: in the best case (met for few transitions), /(i)'''" 
is a good approximation to f^^' and rigidly shifts the whole spectra conserving the shape (and therefore the bound 



Figure 3, Bruneval et al., J. Chem. Phys. 
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FIG. 3: Diagonal of the static inverse RPA dielectric function of bulk silicon. The vertical line denotes the border of the 
Brillouin zone. 

exciton); otherwise, /(i)-''" is numerically instable and gives rise to scattered spectra |30j. Therefore, in the following 
calculations, we always use the non-linearized version ]^^^ (instead of /(i)^''"). For the same reason, it is not astonishing 
that Kim and Gorling have found a "collapse" of the silicon absorption spectrum |lj| . The authors have solved the 
problem by cutting off the long-range (small-g) part of the Coulomb interaction. In fact, this procedure leads to 
drastic changes in the spectrum: we have repeated their calculation by introducing the same cutoff in /'•^^ Now, 
instead of the strongly bound exciton we find a lineshape in good agreement with experiment, as can be seen by the 
dot-dashed curve in Fig.|3 (/'•^-' has not been modified; therefore the spectrum stays in the HF energy region.) In the 
same way, the reduction of the long-range Coulomb interaction in Z'^' translates into a shift of the spectrum towards 
the experimental position. 

This result may seem rather ad hoc. However, it can be understood and used to improve the approach, as we will 
discuss in the following section. 



CORRELATION CONTRIBUTIONS 



Fig. 131 shows the diagonal of the static inverse dielectric matrix eQ^Q(q) of bulk silicon as a function of |q + G|, 
calculated in the Random Phase Approximation (RPA) . The vertical line denotes the border of the Brillouin zone, 
up to which Kim and Gorling 1^ have chosen to set the Coulomb interaction to zero. The step function that one 
obtains in this way can be seen as a first reasonable approximation for the full screening curve. In other words, the 
modified Coulomb interaction is an approximation to the screened Coulomb interaction W: it compensates for the 
lack of correlation. Hence, the impressingly good result of Kim and Gorling and in Fig. |21 can be explained: the 
new, screened p"^' is just an approximation to the electron-hole attraction term derived in Refs. 16, 17 from the 



Figure 4, Bruneval et al., J. Chem. Phys. 
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FIG. 4: Imaginary part of the macroscopic dielectric function of bulk silicon. Continuous curve: non-linearized TD- 
(EXX+cLDA) result; stars: non-linearized TD-EXX result; dot-dashed curve: using /"'•''" and the modified Coulomb in- 
teraction v/6. 



Bethe-Salpeter equation, which it reproduces in the same way as the unscreened version reproduces TD-HF. The 
same applies in principle to f^^'. 

It should be pointed out that the good results obtained from the more sophisticated approaches |l6lll7| like p/G [l5j 
rely in practice on a number of approximations that are commonly made in the Bethe-Salpeter approach from which 
they are derived. In particular, QP eigenvalues are calculated within the GW approximation (including dynamical 
effects), whereas W for the electron-hole screening is taken static. Although these are much less crude approximations 
than the cutoff used above for TD-EXX, the search for a perfectly rigorous, but still efficiently working approach is 
not yet completed. 

At preserit, however, one may be with no doubts satisfied with the precision and reliability of the screened approaches 
|l6Lll7l ll8lll9 | . Nevertheless, it is interesting to investigate the role of correlation further, since the cutoff approach 
of Kim and Gorling gives precious hints: the screening of the long-range (small-g) part of v is seen to have drastic 
effects. This is consistent with other studies of the long-range contribution of the exchange-correlation kernel in bulk 
materials (see e.g. Refs. l8l[l32l) . 

Without considering these findings, one might hope to introduce correlation in another, more standard way, namely 
by adding LDA correlation to EXX as it is quite frequently done for the ground state potential _12J. Fig. 0] shows 
the result of a TD-(EXX-fcLDA) calculation (continuous curve) as compared to the TD-EXX result (stars): the 
LDA correlation does not give rise to any visible changes. (Note that in both cases /(i)-''" has been replaced by the 
corresponding f'-^\ for clarity. In other words, only the effect of correlation in f^^^^ is tested.) 

The adiabatic LDA kernel is in fact short ranged, and cannot suppress the overbound exciton stemming from long- 
range contributions. Instead, as can be seen from the cutoff result, any model e^^ that reasonably screens the long 
range contributions can do the job; for illustration, we also show in Fig. 4 the result obtained replacing v by v/{eo) in 



j(2),im^ where (eo) is taken to be an average dielectric constant of 6 for silicon. In spite of the extreme simplicity of 
the treatment of correlation, the result is again satisfactory (note that, as in Fig.^l the unscreened /^^^ is used; hence, 
the spectrum results too high in energy) . This leaves the hope that, starting from some screened version of TD-EXX 
and along the lines of Refs. [la, llTl Il8l llJA \33L it is possible to find approximations to TDDFT, less rigorous than 
TD-EXX, but exempt from its severe shortcomings concerning the description of bulk materials, and that the method 
is still numerically advantageous with respect to the solution of the Bethe-Salpeter equation. Hybrid functionals like 
the ones discussed in Ref. |3J may be seen as a possibility in this context. 

CONCLUSIONS 

A time-dependent OEP procedure constructs TD-DFT kernels that yield the same time-dependent density as a 
given approximation to the self-energy, via a time-dependent Sham-Schliiter equation. Apart from a linearization, 
this relation is exact. It is therefore not surprising that a careful time-dependent OEP calculation reproduces the 
solution of the Bethe-Salpeter equation, within the corresponding approximation. The present work verified this 
agreement for the case of the Hartree-Fock approximation to the self-energy. In particular, the non-linearized and the 
usual linearized TD-EXX are shown to reproduce the TD-HF calculation and consequently, fail crudely to describe 
absorption spectra of semiconductors, because the electron-hole interaction is largely overestimated there. We show 
that TD-EXX gives rise to a huge bound exciton for bulk silicon. The linearization does not cure this shortcoming. 
The similarity between TD-HF and TD-EXX can also be noticed in the evaluation of vertical excitation energies of 
finite systems J35]. 

Whatever the approach used, Bethe-Salpeter equation, p/G method, or TD-DFT, the inclusion of the screening of 
the exchange operator is evidenced as crucial to give a proper account for the electron-hole interaction. In particular, 
the long-range components of the exchange kernel have to be reduced. Furthermore, we demonstrate that the details 
of the screening do not matter much: two very different and crude models (a cutoff as used by Kim and Gorling or a 
uniformly reduced Coulomb interaction w/6) allow us to equally produce realistic absorption spectra of silicon. 

The hope of describing optical absorption spectra of semiconductors within TD-DFT is well founded. However, in 
the framework of unscreened methods, like TD-EXX is, there is no chance to get something else than the disastrous 
TD-HF results: one has to go beyond TD-EXX. A rigorous OEP method based on "exact screened exchange" may 
do the job. Fortunately, the inclusion of screening within rather simple approximations, e.g. using an empirically 
screened Coulomb interaction, seems to be already sufficient. 

We thank Andrea Cucca for his help. We are grateful for discussions with R. Del Sole and A. Gorling, and for 
computer time from IDRIS (project 544). This work has been supported by the EU's 6*'' Framework Programme 
through the NANOQUANTA Network of Excellence (NMP4-CT-2004-500198). 

DERIVATION OF THE TD-DFT KERNELS FROM THE LINEARIZED TD-SHAM SCHLUTER 

EQUATION 

The present appendix provides the derivation of the linearized TD-EXX kernels from the linearized TD-Sham- 
Schliiter equation. The linearized TD-Sham-Schliiter equation reads 

/dS Gexx(1,3)«exx(3)Gexx(3,1) = f d3d4 Gexx(1,3)E,(3,4)Gexx(4, 1). (14) 

When Eq. 114(1 is differentiated with respect to the TD density p{2), we get 

J d3 Gexx(1, 3)^^|^^Gexx(3, I) = J dUA Gexx(1, 3)^|^^Gexx(4, 1) 

+ j dUA '^^y^^^l'^^ [S.(3, 4) - 5{i, 4)«Exx(3)] Gexx(4, 1) 

+ j didA Gexx(1, 3) p,(3, 4) - ,5(3, 4)«exx(3)] '^^7^2!' ^^ ^ ^^^^ 

The linearized exchange operator is simply Sa;(l,2) = iGEXx(l, 2)11(1, 2). Therefore, the only quantity needed to 
carry on the derivation is the derivative of Gexx with respect to p. It can be evaluated along the following lines. 
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using standard functional analysis relations, and introducing the total KS potential within EXX i;kS: 
<5Gexx(1,2) ^ f ^^ 5Gexx{1,2) SvKsi^) 



SpiS) J 6vKsi^) Sp{3) 



d4d5rf6GExx(l,5)GExx(6,2)- 

/ 



dA Gexx(1, 4)Gexx(4, 2)x^''''"'(4, 3), (16) 



where the last line was obtained from the Dyson equation Gg^x = '^o^ ~ "ks (Gq standing for the free-electron 
Green's function) and from the definition y^^^^ = Svks/^P- 

Finally, by inserting Eq. Hlf)|l into Eq. H15|) . by multiplying by Xq'^-^(2, 5), and integrating over the variable 2, we 
obtain the central equation for the linearized TD-EXX kernel /^xx, im _ Svexx/^P- 

J d2di Xo''''(l, 3)/^'''''''"(3, 2)x|;'^^(2, 5) = y" didA Gexx(1, 3)Gexx(4, 1)i'(3, 4)Gexx(3, 5)Gexx(5, 4) 
- zGexx(1, 5) j dMA Gexx(5, 3) [S,(3, 4) - <5(3, 4)t;Exx(3)] Gexx(4, 1) 

- iGexx(5, 1) / dUA Gexx(1, 3) [I],(3, 4) - ,5(3, 4)z;exx(3)] Gexx(4, 5), (17) 



with Xo^^ = -^GexxGexx- The kernel f^^^^ "" can be split into two pieces /(i)>i™ and Z^^^'"", the definitions of 
which stand respectively in Eqs. IjlSfl and H12(l . This partition is natural when looking at the analytical form of the 
terms. It is further physically-driven, since the term /(2)'''" accounts for electron-hole interaction and the term /f^)^''" 
for the quasiparticle shift. I1!tLI28I| 

LINK TO TD-EXX 

This appendix shows that the kernel obtained in the previous appendix is precisely the EXX kernel of Kim and 
Gorling "29). For simphfication, let us name T^^")-''" the first term of Xo^^/^^^'""Xo^^' T^^b).\\n ^j^g second one, and 

7^(2) Jin „ EXX ^(2) Jin EXX 
AO J AO 

We are about to introduce the expression of Gexx in the previous terms in order to recover all the sixteen terms 
of Kim and Gorling's kernel. The time-ordered EXX Green's function Gexx in frequency domain is 



cj)^{Yl)(j)*{v2) 



GExx(ri,r2,c^) = >^ "'^ ."' \^ ,, , (18) 






w - e. 



z77(2/,-l)' 



where (pi and e^ are the EXX KS wavefunctions and energies for index i (that contains also the k point information) . 
/,; is 1 for occupied states and for empty states. 

Evaluation of T^^)"" 

Let us first proceed with the electron-hole interaction term T^^^'^'". It reads, after Fourier transform to frequency 
domain, 



r(2)-'-(ri,r5,c.) = ^^0,;(ri),/);(ri) j du,^- \ ^ 



ijkl 

X {tk\v\jl)rkir5)M^5) fdu^2^— ^7 r, (19) 

J (w + W2 -efc)(w2 -£/) 

as the products in time space become convolutions of frequencies. The factor 2 accounts for spin degeneracy. The 
usual Coulomb integrals 

{ik\v\jl)^ j dv^dv2<l^*{y^)<j,k{y^)—^<P,{v2Wi{v2) (20) 

J n ~ 1^2 
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have been introduced and the Hrj factors in the denominators are still present, but not explicitly written (they are 
unchanged with respect to the definition of Gexx)- 

The frequency integrals are now calculated by virtue of the residue theorem on a path that encloses cither the 
upper half-plane, or the lower half-plane. Contributions with all poles in the same half-plane vanish. Consequently, 
the frequency integrals are 

J ^ n ■ Jj Ji (OT) 

^[uj + uji- Ei - i?7(2/j - l)][wi - Ej - ir]{2fj - 1)] w - {e, - e^) + irjifj - fi) ' 

The t(2)'"" term finally reads 

r^ ^' (ri,r5,w) = -2 > (/j -/i) ^ , . .„ -p-{ik\v\jl){ii- }k) -, \ , ■ if FT- (22) 

This expression for T^^^^''" is equal to the H^ and Hj^ terms of Ref. 123, except that the convergence factors irj 
are of opposite sign for antiresonant terms. In fact, the present derivation, starting from time-ordered Green's 
functions, yields time-ordered quantities, whereas Kim and Gorling's derivation considers causal quantities. When 
used adequately, this difference is not relevant in practical applications. 

Evaluation of T'^"''"" and T^"'"" 

Let us now turn to the contribution ^(i'')^^'" to the linearized TD-EXX kernel. 'Sx is a static approximation for the 
self-energy, hence in the frequency domain, T^^^"^'''" reads 

2tt ^ j uj + uji-e,- ir?(2/i - 1) 

X —7777 TriJl^x -VEXxW -7:77 TT- i^-^j 

uji - e-j - iri[2jj - 1) uji- Ek- ir]{^Jk - 1) 

Performing the integration on ui thanks to the residue theorem gives a vanishing contribution if the i, j, k states 
are all occupied or all empty. There are 6 non-vanishing terms corresponding to the other cases. We will exemplify 
three of them in the following. The three remaining ones are analogous. 

If i and j are occupied and k is empty, let us close the path of integration in the lower half-plane. The enclosed 
poles are the ek — irj that yield the residues: 

-Zm- —-, y (24) 

[uj- [ei - ek) -iri\[ek - ej\ 

If i and k are occupied and j is empty, closing the path analogously in the lower half-plane retains the poles ej — irj 
that give the residues: 

_2m /»(! - fj)fk ^25) 

If J et k are occupied and i is empty, this retains poles located at e^ — a; — ivy with residues: 

_2_-v (1 ~ JijJjJk ,„„•. 

[uj - (e^ - ej) - ii]][uj - (e^ - ek) - irj] ' 

The three other terms correspond to the case with 2 empty states and 1 occupied. The path of integration will be 
closed in the upper half-plane, in order to retain only the poles from the occupied states. 

j-Cia),!"! g^ally gives rise to 6 terms. T'^^''^''™ will also account for 6 analogous terms. The sum of T^^^^'^™ and 
rp(ib),hn^ if written explicitly, is exactly the terms H^+H\ of Kim and Gorling, except once again that the convergence 
factors irj are opposite for antiresonant transitions. 

This appendix showed that the linearized Sham-Schliiter equation indeed yields the same TD-DFT kernel, as the 
one obtained by Kim and Gorling in Ref. |29|. 
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